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The spin-Peierls transition in the ground state of a quasi-one-dimensional spin-i Heisenberg model 
coupled to adiabatic bond phonons is studied using a quantum Monte Carlo (QMC) method. The 
transition from a gapless Neel state to a spin-gapped Peierls state is explored in the parameter space 
spanned by the spatial anisotropy and the strength of spin-lattice coupling. It is found that for any 
finite inter-chain coupling, the transition to a dimerized Peierls ground state occurs only when the 
spin-lattice coupling exceeds a finite, non-zero critical value. This is in contrast to the pure ID 
model, where adiabatic phonons lead to a dimerized ground state for any non-zero spin-phonon 
coupling. The phase diagram in the above parameter space is mapped out. No evidence is found 
for a region of co-existing long range magnetic order and dimerization. The vanishing of the Neel 
order occurs simultaneously with the setting in of the dimerization. 



PACS: 75.40.Gb, 75.40.Mg, 75.10.Jm, 75.30.Ds 
I. INTRODUCTION 

When coupled to an elastic lattice, a spin-i Heisen- 
berg chain is unstable towards a dimerization of the lat- 
tice. The cost in elastic energy due to a distortion 5 of 
the lattice 5"^) is compensated for by the gain in the 
magnetic energy (~ 5'^^'^), and the ground state is sta- 
bilized for a lattice with non-zero dimerization^^^ and a 
finite spin gap. Such a transition is known as the spin- 
Peierls (SP) transition - in analogy with the conventional 
Peierls transition in a one-dimensional (ID) metal. The 
SP transition has been extensively studied and is believed 
to be fairly well understood in 1D.^~^ In the adiabatic 
limit, any arbitrarily small coupling to an elastic lattice 
leads to a dimerized ground state. On the other hand, 
quantum lattice fluctuations destroy the dimerization for 
small spin-phonon couplings and/or large bare phonon 
frequencies.^'^ The transition to a dimerized (Peierls) 
state occurs only when the spin-phonon coupling ex- 
ceeds a finite, non-zero critical value (that depends on 
the bare phonon frequency). This is true for both the 
spin-i Heisenberg and XY models with spin-phonon cou- 
pling. The discovery^*^ of the quasi-lD inorganic spin- 
Peierls compound CuGeOa has led to a resurgence in the 
study of spin-Peierls(SP) transition in low-dimensional 
spin models. The properties of CuGeOs have been widely 
studied within the framework of a ID spin-i Heisenberg 
model coupled to phonons, with an additional frustrated 
next-nearest-neighbor interaction.^^ 

Fewer studies exist for the spin-Peierls transition in 
two dimensions (2D) and the mechanism of the transi- 
tion is not as well established as in ID. Although the 
study of the Peierls transition in 2D or quasi-lD sys- 
tems presents a technically much more difficult challenge, 
it also allows for a more realistic modeling of real ma- 
terials - for CuGeOs, the strength of inter-chain cou- 
pling is estimated to be J±/J ~ 0.1 experimentally^^, 
and its importance in any realistic modeling of CuGeOa 



has been widely stressed. From a purely theoretical 
standpoint, the study of spin-Peierls transition is inter- 
esting on its own. Unlike in ID, the ground state of the 
2D (Heisenberg) system, in the absence of spin-phonon 
coupling, has long range antiferromagnetic (Neel) order. 
It is generally believed that even for adiabatic phonons, 
the spin-phonon coupling has to exceed some non-zero 
critical value for the ground state to develop a dimer- 
ized pattern with a spin gap. What is the nature of the 
transition? Is there a region in the phase space where 
the ground state has co-existing dimerization and long 
range antiferromagnetic order? Furthermore, the exis- 
tence of several different possible dimerization pattern in 
2D means that, in principle, different dimerization pat- 
terns can be stabilized for different values of the param- 
eters. Finally, unlike in ID, the Peierls phase in 2D ex- 
tends to finite temperatures. 

The spatially isotropic 2D spin-i Heisenberg model 
with static dimerization patterns have been studied by 
several authors. By comparing the ground state en- 
ergies for different patterns, these authors have tried 
to predict the energetically most favored dimerization 
pattern. Using the same strategy, the 2D tight bind- 
ing model with bond-distortions^'^ and the 2D Peierls- 
Hubbard model'^^^'^^ have also been studied. Since in the 
limit of large on-site repulsion, [/, the Hubbard model 
at half-filling reduces to the Heisenberg model, results 
from the Peierls-Hubbard model (in the limit of large 
U) are expected to be applicable to the present discus- 
sion. Unfortunately, there is no consensus among the 
different studies as to the nature of the dimerization 
pattern in the ground state. For the Peierls-Hubbard 
model at half-filling. Tang and Hirsch^^ find a plaquette- 
like distortion to be energetically favored in the limit 
of large U . On the other hand, Mazumdar^"^ have ar- 
gued that the minimum energy ground state has a "stair- 
like" dimerization pattern - corresponding to a wave vec- 
tor Q = (tt, 7r). Zhang and Prelovsek^^ agree with a 
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dimcrization pattern with Q = (tTjTt), but conclude that 
the ground state has dimcrization only along one of the 
axes (staggered dimerized chains - similar to the pat- 
tern considered here). For the Heisenberg model with 
static dimerization, Al-Omari^^ has concluded that the 
ground state energy is minimized for a plaquette-like dis- 
tortion of the lattice which agrees with the conclusion of 
Tang and Hirsch. On the other hand, Sirker et.al?^ find 
that linear spin wave theory (LSWT) predicts a stairlikc 
dimerization pattern to be most favored, in agreement 
with Mazumdar. Using LSWT, Sirker et.al. also find fi- 
nite regions in the parameter space with co-existing long 
range magnetic order and non-zero dimerization. How- 
ever, as pointed out by the authors, results obtained from 
LSWT are not reliable at large values of dimerization. 
The effects of inter-chain coupling was considered early 
on by Inagaki and Fukuyama^^'^^ who studied a quasi- 
ID system of coupled spin-i Heisenberg chains with a 
fixed dimerization pattern corresponding to a wave vec- 
tor Q = (7r,0). By treating the inter-chain coupling in 
a mean-field theory, they were able map out the ground 
state phase diagram^^ and study the finite temperature 
transition-'^^ . Later, Katoh and Imada-*^^ studied in detail 
the nature of the transition. More recently, the effects of 
impurities have been studied in this model (once again 
with a mean-field treatment of the inter-chain coupling) 
which revealed a region of co-existing Peierls and anti- 
ferromagnetic orders. ^^""^^ In addition to this, the quasi- 
ID XY model with a Q = (7r,7r) dimerization pattern 
has recently been studied^^'^^ by an extension of the 
Jordan- Wigner transformation in 2D.^'^~"^''' The effects of 
quantum phonons in the isotropic 2D model have also 
been studied,^^ where the authors find that there is no 
evidence of a transition to the Peierls state for a wide 
range of values of the bare phonon frequency and the 
spin-phonon coupling. This is consistent with previous 
finding^'' for the same model that the spin wave spec- 
trum along the Brillouin zone boundary is qualitatively 
similar to that for the pure Heisenberg model. 

The present work aims to investigate the SP transition 
in a spin-i quasi- ID Heisenberg antiferromagnet coupled 
to static Q = (tt, it) bond phonons with varying inter- 
chain coupling. The elastic energy due to the lattice 
distortion is considered as well as the change in mag- 
netic energy. The nature of the ground state is explored 
for different values of the inter-chain coupling and the 
elastic constant associated with the bond distortions and 
the ground state phase diagram in this phase space is 
mapped out. 

The rest of the paper is organized as follows. In Sec- 
tion II, the model Hamiltonian and the Stochastic Series 
Expansion (SSE) QMC method used to study it are in- 
troduced. The results of the simulations are presented in 
Section III. Section IV concludes with a summary of the 
results. 
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FIG. 1. The spin part of the ground state energy as a func- 
tion of the system size for four values of 6, and a=0.25. 



II. MODEL AND SIMULATION TECHNIQUES 

A quantum Monte Carlo (QMC) method has been used 
to study a quasi-lD Heisenberg model with spin-phonon 
coupling. The model is given by the Hamiltonian 



+ J± ^ Sij ■ Sjj+i , (1) 

hi 

where J± is the inter-chain coupling, A is the strength 
of the spin-phonon coupling (restricted to be only along 
the chains), Uj.j's arc the distortions of the bond lengths 
and K is the elastic constant for the distortions. Follow- 
ing the experimentally observed^^ dimerization pattern 
in CuGeOa, the bond length distortions are chosen to be 
of the form 

This amounts to choosing a fixed dimerization pattern 
corresponding to the wave- vector Q = (7r,7r). As in 
previous works, the bond distortions can be rescaled by 
the spin-phonon coupling strength, A. This reduces the 
Hamiltonian to 
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FIG. 2. The spin part and the total ground state energy 
as a function of the bond distortion for 3 different values of 
the elastic constant, and fixed a=0.25. The system size is 
iV=64xl6. 

H = ^(1 + (-1)^+^(5)8, J • Si+i,, + Ndy2C 
hi 

+aJ2Sij-Sij+i, (2) 

where C = ^j^, a = Jj_/ J and N is the size of the lattice. 
The static approximation for the displacements makes 
the computational task easier. The following strategy is 
adopted. The simulations arc carried out for the spin 
variables for several different {a, 6} parameter sets. This 
produces the spin energy of the system as a function of 6 
for a fixed a. Next the elastic energy with a particular C 
is added and the total energy is minimized to obtain the 
value of the ground state distortion for the given set of 
parameters {a, (}. This is repeated for different sets of 
{a, (} to obtain the ground state phase diagram in the 
parameter space spanned by a and (. 

The stochastic series expansion (SSE) QMC method 
has been used to study the present model. The SSE 
method^^"^^ is a finitc-tcmpcraturc quantum Monte 
Carlo method based on importance sampling of the di- 
agonal elements of the Taylor expansion of e~^^ , where 
(3 is the inverse temperature /? = J/T. Ground state ex- 
pectation values can be obtained using sufficiently large 
values of /3, and there are no approximations beyond sta- 
tistical errors. Using the recently developed "operator 
loop update" it is possible to explore the phase space 
of the Hamiltonian (2) in an efficient manner. 

III. RESULTS 

For a system of weakly coupled Heisenberg chains, it 
was shown recently that the estimates for various ob- 
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FIG. 3. The equilibrium ground state bond distortion as a 
function of the spin-lattice coupling for different values of the 
inter-chain coupling. 



servablcs for a spatially anisotropic system depend non- 
monotonically on the system size for square {L^ = Ly) 
geometry.*^ One has to go to rectangular (L^ i= Ly) ge- 
ometries to obtain monotonic behavior of the numeri- 
cal results for extrapolating to the thermodynamic limit. 
Similar effects are expected in the present model for 
a <C 1. Hence rectangular lattices with the aspect ratio 
Lx = 4Lyhave been studied, with = 16 — 64. An in- 
verse temperature of /3 = SL^ was found to be sufficient 
for the observables to have converged to their ground 
state values. Nine different values of the inter-chain cou- 
pling wore considered, a = 0.01, . . . , 1.00. 

Fig.l shows the finite-size dependence of the QMC 
data for the spin part of the ground state energy for four 
different values of the bond length distortion, 5. The 
values of 5 chosen correspond to a Neel ground state 
((5=0.04, 0.10), a Peierls ground state (5=0.20), and a 
groimd state close to the critical 5 for Nccl to spin-Pcicrls 
transition in the pure spin system (the phase transition 
in the spin part of the Hamiltonian (2) is discussed at 
the end of the current section). The leading order finite- 
size correction to the ground state energy is seen to be 
^ l/iiJ for small values of i5 corresponding to a Neel- 
ordered ground state. This is similar to that observed for 
the pure 2D Heisenberg model. *° On the other hand, for a 
spin-gapped groimd state ((5=0.20), the energy scales ex- 
ponentially with system-size. Close to the critical point, 
extrapolation to the thermodynamic limit becomes diffi- 
cult due to cross-over effects. Instead, the data from the 
largest system size studied have been used to map out 
the phase diagram. Fortunately, the data for the largest 
system sizes studied are found to be well converged the 
fractional difference in the energy for the two largest sys- 
tem sizes studied is ~ 10~^. This observed convergence 
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allows for reliable estimation of ground state properties 
in the thermodynamic limit based on the data from the 
largest system size any finite-size effects on such esti- 
mates are expected to be small. 

The strategy implemented to extract the ground state 
bond distortion is qualitatively demonstrated in fig. 2. 
The total ground state energy is obtained by adding the 
elastic energy contribution to the spin part of the energy 
obtained from the simulations. The plot shows the spin 
part of the energy, as well as the total ground state en- 
ergy as a function of the bond length distortion, 5, for 
three different values of the clastic energy constant, C, 
and a fixed value of the inter-chain coupling (a^O.25). 
For large C, the total energy is minimized for a non-zero 
value of the bond length distortion, S, whereas for small 
C, a uniform ground state with 5=0 is energetically fa- 
vored. The behavior of the total energy near the critical ( 
is also shown. The ground state distortion in bond length 
is obtained by numerically differentiating the total en- 

ergy data and solving for \sgs = 0. In principle, one 

can also fit a polynomial to the QMC data to get E'*^™ (S) 
and add to it the elastic energy term to get E*°*(J). Then 
the ground state distortion can be obtained as a contin- 
uous function of S by solving g^. \g93 = 0. However, in 
practice, the earlier method is found to be more reliable 
because of the uncertainty in the order of the polynomial 
fit. 

Fig. 3 shows the equilibrium distortion in the bond 
lengths in the ground state of the system as a fmic- 
tion of the elastic energy parameter, at fixed values 
of a, obtained as described above. For small values of (, 
the tendency towards dimcrization is suppressed, and a 
uniform (Neel ordered) ground state with uniform bond 
lengths is stabilized. As ( is increased above a critical 
value, Cc. there is a transition to a ground state with 
a finite, non-zero dimerization. For ^ > the equi- 
librium ground state distortion increases monotonically 
with The critical value Cc required for a transition to a 
dimerized ground state depends on the value of the inter- 
chain coupling, a and increases with increasing a. The 
data indicate that the transition is a discontinuous (first 
order) quantum phase transition for any finite value of 
inter- chain coupling. 

The results from fig. 3 are summarized in fig. 4, which 
shows the phase diagram for the system in the phase 
space spanned by the parameters ( and a. For small C 
and/or large a, the ground state of the system is Neel or- 
dered with zero bond distortion, while for large C and/or 
small a, the ground state is dimerized with a finite spin 
gap. The phase boundary is roughly linear over the range 
of values of a studied. However, it is well known that in 
the ID limit (a = 0), the critical elastic constant should 
vanish (Cc = 0). Unfortunately, the numerical data be- 
comes too noisy in this limit for an accurate estimate 
of the critical coupling. For the quasi- ID XY model, it 
was shown that Cc goes as — 1/lna for a — > 0. On the 
other hand, results for a system of coupled (undimerized) 
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FIG. 4. The ground-state phase diagram in the parameter 
space of the inter-chain coupling and the spin-lattice coupling 
strength. 

Heisenberg chains show that the magnetization vanishes 
as a power law in a with a logarithmic correction.^^ A 
similar behavior is expected for the critical coupling in 
the present model. 

This brings up the interesting possibility of having a 
region in the (C, a) parameter space where the ground 
state has co-existing Neel order and non-zero dimeriza- 
tion. Such co-existence has been shown to exist in the 
presence of doping. ^^"■'^ For this purpose, the static spin 
susceptibility defined as 

^(q) = ^E^''''^"''"''^^ dr{S^{T)Stm. (3) 

have been studied for the spin part of the Hamiltonian 
(2) (without the elastic energy term). 

hi 

hi 

If the ground state has long range antiferromagnetic or- 
der, the staggered (Q = (tt, tt)) susceptibility scaled by 
the system size {x{tt,'jt)/N), for a finite system will in- 
crease with increasing system size, diverging in the ther- 
modynamic limit. On the other hand, if the ground state 
has a finite spin gap, x(7r, 7r)/A'' will vanish in the limit 
of infinite system size. This qualitative criterion can be 
expressed in a more quantitative manner by noting that 
the ground state of the above Hamiltonian undergoes a 
continuous transition from a Neel ordered state with long 
range antiferromagnetic order to a spin-gapped, dimer- 
ized phase with no long range magnetic order as 6 is in- 
creased beyond a finite, non-zero critical value, 6*, that 
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depends on the inter-chain coupling a. The transition 
belongs to the universality class of the 3D Heisenberg 
model'*'^. Finite-size scaling*^ predicts that for such a 
transition, the finite-size susceptibility at the critical 5 
scales with the system size as 

X(L.) ~ Ll-\ 

for a rectangular lattice of dimension N = L^xLy. This 
implies that on a plot of x(7r, 7r)/i^^'', the curves for dif- 
ferent system sizes will intersect at the critical S. The 
value of the critical exponent rj is known to a high degree 
of accuracy {r] « O.OS?)."*^ Fig. 5 shows x(7r, 7r)/L^~'' as 
a function of S for a fixed value of a = 0.25 for several 
different system sizes. For small 6, the scaled susceptibil- 
ity increases with increasing system size, indicating the 
presence of long range magnetic order. For larger values 
of (5, the scaled susceptibility goes to zero with increas- 
ing system size, signaling the opening up of a spin gap. 
From the data, the critical value of S is estimated to be 
S* « 0.16. This value of the bond length distortion is 
less than the jump in 6 at the transition point, Cc. This 
means that for C > Cc) the ground state is dimerized with 
§gs J* g^jj(j jjj^g jj^Q jQjjg range magnetic order. On the 
other hand, for C < Cc, the a uniform ground state is en- 
ergetically favored that has zero dimerization {S^^ = 0), 
and long range magnetic order (the Neel state). For no 
values of C is a ground state with < S^^ < S* ever fa- 
vored energetically. Hence the transition to the dimerized 
phase is accompanied by the disappearance of magnetic 
order and there is no region of co-existing dimerization 
and magnetic order. This is true for the present model. 
Whether it is possible to have ground states with co- 
existing magnetic order and dimerization in other mod- 
els (with different dimerization patterns) remains to be 
seen. 



IV. SUMMARY 

A quantum Monte Carlo method has been used to 

study a quasi-lD Heisenberg model coupled to static 
bond phonons. Motivated by experimental observa- 
tions in the inorganic quasi-lD spin-Peierls compound 
CuGeOa, the bond distortions are restricted to be only 
along the chains, and a single dimerization pattern, cor- 
responding to the wave vector Q — (tt, tt), is considered. 
It is found that unlike the pure ID case, for any non-zero 
inter-chain coupling a, the transition to a dimerized spin- 
Pcierls ground state occurs only when the spin-lattice 
coupling strength, (, exceeds a finite, non-zero critical 
value, Co even for the static phonons considered here. 
For C < Cc, the ground state has long range Nccl order 
and zero spin gap, whereas for C > Cc, the ground state 
develops a finite dimerization accompanied by the open- 
ing up of a spin gap. The transition is found to be a dis- 
continuous (first order) quantum phase transition. The 
value of the critical coupling depends on the strength of 
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FIG. 5. The ground-state static staggaod magnetic sus- 
ceptibility as a function of bond distortion for a fixed value 
of the inter-chain coupling, a = 0.25. 

the inter-chain coupling, increasing monotonically with 
a. The phase diagram in the parameter space of a and 
C is mapped out. Finally it is found that in the present 
model, the transition to the dimerized Pcicrls state is ac- 
companied by the vanishing of magnetic order, and that 
there is no region of co-existing magnetic order and non- 
zero dimerization. 
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